By counting the number of unique species are found in each of the half degree cells we can then create a single global raster of number of species.
## quartz_off_screen
## 2
## quartz_off_screen
## 2
| class | AquaMaps | IUCN | Both |
|---|---|---|---|
| Actinopterygii | 11156 | 5540 | 1460 |
| Anthozoa | 788 | 577 | 394 |
| Ascidiacea | 536 | NA | NA |
| Asteroidea | 161 | NA | NA |
| Aves | 188 | NA | NA |
| Bivalvia | 833 | NA | NA |
| Cephalopoda | 246 | NA | NA |
| Demospongiae | 344 | NA | NA |
| Echinoidea | 113 | NA | NA |
| Elasmobranchii | 809 | NA | NA |
| Gastropoda | 2270 | 740 | 125 |
| Holothuroidea | 78 | 396 | 44 |
| Hydrozoa | 234 | 10 | 2 |
| Magnoliopsida | 5 | 64 | 4 |
| Malacostraca | 2897 | 302 | 161 |
| Mammalia | 118 | 156 | 109 |
| Maxillopoda | 175 | NA | NA |
| Ophiuroidea | 150 | NA | NA |
| Polychaeta | 204 | NA | NA |
| Pycnogonida | 466 | NA | NA |
| Reptilia | 32 | 210 | 27 |
| Scaphopoda | 223 | NA | NA |
| Liliopsida | NA | 89 | NA |
| Polypodiopsida | NA | 3 | NA |
This looks at the full 2014 Aquamaps dataset… should it be limited to species with valid categories?
## Joining by: c("sciname", "iucn_sid", "am_sid", "occurcells")